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Abstract. 

A distribution of electromagnetic fields presents a statistical assembly of a particular 
type, which is at scale h a quantum statistical assembly itself and has also been instrumental 
to a concrete demonstration of the basic probability assumption of quantum mechanics. Of 
specific concern in this discussion is an extensive wave train of radiation fields, described by a 
total wave function ip, which are continuously (re)emitted and (re)absorbed by an oscillatory 
(point) charge of a zero rest mass and yet a finite dynamical mass, with the waves and charge 
together making up an extensive undulatory IED particle. The IED particle will as any real 
particle be subject to interactions with the environmental fields and particles, hence to 
excitations, and therefore will explore all possible states over time; at scale h the states are 
discrete. On the basis of the principles of statistics and statistical mechanics combined with 
first principles solutions for the IED particle, we derive for the specific statistical system of 
the IED particle the probability functions in position space, of a form IV^ 2 ! and in dynamical- 
variable space. 



1. Introduction 

It has been a long established fact in experiment and formal theory that matter particles 
manifest both the characteristics of corpuscle and wave, termed as matter waves (<?■), and 
are at the scale of Planck constant h dominated by quantum mechanics which formally 
is formulated on the basis of the wave equations of E Schrodinger, W Heisenberg and P 
Dirac[l]. The interpretation of the quantum mechanical wave function ^ is based on the 
probability assumption of M Born [2] that the state of a particle, assumed corpuscular, is 
completely specified by its ^(x, f), with ^(x, t) being generally complex, and \\P(x, t)\ 2 represents 
the statistical probability of finding the particle at position x. This probability has a key 
feature distinct from that of classical systems: it is additionally constrained by the Heisenberg 
uncertainty relations for the conjugate dynamical variables [3]. In phase space, this probability 
is associated with a smallest volume, equal to the fundamental constant h of M Planck raised to 
power fN for a N particle system of / degrees of freedom, which is accessible to a microscopic 
state. 

Quantum mechanics has proven successful in predicting a broad range of quantum 
mechanical properties of particles, especially properties involving the quantization of dynamical 
variables, for which the quantum theoretical predictions are directly comparable with the overall 
vast quantity of experimental data accumulated over a century or so to this date. However, it 
is up to the present not satisfactorily understood that what is waving with a matter wave, what 
is accordingly statistical with it in position as well as dynamical variable spaces, and what is 
the origin of the Planck constant which defines a smallest volume in phase space? 

Their understandings come especially in question when our purpose is in order that the 
probability, and the fundamental constant h here, can be inferred based on a common minimal 
set of first principles' laws which also govern the remaining basic properties of particles. Further 
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to this, quantum mechanics as interpreted through a statistical corpuscular-particle picture 
meets difficulties in (coherently) accounting for certain basic experimental phenomena, including 
in particular the basic coherent wave phenomena such as diffraction and (self) interference [4]. 
Furthermore, the corpuscular particle with a given rest mass does not furnish a mechanical 
scheme for accounting for the origin of mass, the cause of gravity and the origin of relativistic 
effect, among others (see [4] for a recent fuller review). 

On the whole, we appear to have in hand a comprehensive, rigorous mathematical structure 
for quantum mechanics and yet an incomplete conceptual grasp of the underlying mechanical 
content. Conceptually, we are confronted with various internal clashes. These concern for 
examples the dual corpuscle and wave aspects in a single object, the desired wave form of 
spreading and the corpuscular particle assumption at any one time, and the relationships 
between prediction, measurement and actuality of the state of a particle as a subject of debate 
ever since the foundation of quantum mechanics. The difficulties and clashes, as one would 
identify upon an examination, are the inevitable results of descriptions based on the corpuscular 
particle picture. Although, it can not be overemphasised that, the current description represents 
the best theory which is as accurately as can be made in the given framework. 

From the standpoint that nature constructs and operates itself by mechanical means, 
indefinitely down to as small a scale as it takes to sum up to the lager-scale manifestations 
in question, one may anticipate, as it has already been demonstrated to certain extent in [5a- 
m], that the conceptual clashes will all of a sudden vanish once a realistic model of particle, 
one with a realistic internal structure and dynamical scheme, is constructed. 

Based on overall experimental observations as input information, the experimentally 
established de Broglie wave and relations from the outset, the author developed in [5a-m] 
an internally electrodynamic (IED) particle model, which briefly states that a single- charged 
material particle, like the electron, proton, etc., is composed of (i) an oscillatory point-like 
charge q (as source) with a characteristic frequency Q and zero rest mass, and (ii) the 
electromagnetic waves E, B generated by the oscillating charge. The IED particle thus defined is 
geometrically an extensive object, where the term "particle" stands no longer for a "geometrical 
point". Nevertheless, it maintains the other usual aspects pertinent of a "particle", such that 
it separates one distinct kind of material species from others. 

The IED process is itself built on the experimental facts as reflected through a minimal 
set of established basic, or, first principles laws (see [4] for a recent review). Based on the 
first principles solutions for the IED process, it has become feasible to predict a range of basic 
properties of particle as demonstrated in [5a- m] . 

It has been one basic solution obtained earlier [5a,c,d] that the total energy e of an IED 
particle is in direct proportion to the quadratic of the complex total wave function 2p(x,t), i.e. 
e oc \ip(x,t)\ 2 , where \tp(x,t)\ 2 plays an apparent role of probability in position space. In this 
paper we shall deduce the \ip(x,t)\ 2 as a probability in position space more rigorously, and in 
turn its counterpart in dynamical- variable space, on the basis of the basic principles of statistics 
and statistical mechanics, the first-principles solutions for single and many IED particles, and 
the relevant experimental indications. 

2. Wave mechanics of single IED particle. Probability density in position space 

Consider that an IED particle travelling at a velocity v along the x axis between two non- 
absorbing massive walls at x — and L in the presence of an external potential field y(x). 
Its constituent point-like charge q, as source, is endowed with a total mechanical energy (or 
Hamiltonian) e q at initial time, and will over time oscillate continuously at a characteristic 
frequency i?, and also travelling at velocity v in the x direction. Two opposite travelling, 
Doppler-displaced electromagnetic waves of fields W , B J , which together with the q constitute 
the IED particle, are generated by the charge along each direction, where j = f indicates the 
component wave travelling in the direction parallel with v, and j = $ antiparallcl with v. 
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In regions excluding the charge, assuming that, except for y(x), no other charges and 
currents present, the W , B J fields are subject to the Maxwell's equations 

V • W = 0, V • W = 0, V x B j = + \d t W, VxE 3 = -&B 3 ', (1) 

c' 

where 9 t = J^, and c' 2 = c 2 +Y /m, with c(= y/l/eofio) the velocity of light in a free vacuum 
and m the particle's dynamical mass expressed in (6) below. Setting E^> = E q Lp\ with E q the 
amplitude and ^'s the (real) dimensionless fields, or wave functions, their superposition gives 
the real part of a complex total wave function tp = E/E g , Re[tp] = J2j ^ l^q 



Further restricting for illustration here the waves to be propagated along the x axis only, we 
obtain from the Maxwell's equations a total wave equation = (c 2 + ^ {x) / m)V 2 tp . The 
equation, after routine transformations under conditions of small "V and an insured existence 
of stationary-state solutions [5c] , reduces to 

itp = Htp, i = ihd t , H = f/m + y{x), p={h/i)V. (2) 

If after N v cycles of continuous oscillations the charge had emitted its entire initial energy 
s q without reabsorption, nor with supply of external energy, e q is then converted entirely to the 
energy of total wave tp generated by it, e = y/e^e* — e q . The total wave tp extends as a wave 




train of a mean total length L v = - 

In the presence of reflection walls assumed earlier, the radiation waves </5 t ,</? t , and thus 
the total wave tp, will be reflected back from the walls to the charge, be reabsorbed and re- 
emitted by it, repeatedly. Assume that the particle has reached an equilibrium state in [0,L], 
with L being relatively small. Its total energy (e to t) will ordinarily be carried in part by the 
charge and in part by the wave, by the fractions oi and 02, with 01,02 < 1, 01+02 = 1, and 
£tot — oi£ 9 + 02£. The extensivness of the IED particle is determined by its wave components, 
the total tp in question here, and is of our primary concern below, tp will be based on in the 
discussions below. 

Our primary concern here is " probability" . We thus start by directly writing down the total 
electromagnetic energy density, in terms of its total field E(x, t) — E q tp(x, t), in the usual form 
eo(x,t) = eo\E(x, t)\ 2 . e(t) = J Q v eo(x,t)dx = ^f- J £o(x,t)dx gives the total wave energy at 
time t. It is here that the complex form of tp finds its clear physical significance: The complex 
tp ensures that e is a mechanical energy, or, Hamiltonian oc \E\ 2 = (Re[E]) 2 + (lm[E]) 2 7 which 
includes both an inertial component oc (Rc[_E]) 2 carrying the kinetic energy of the field E, and 
an elastic component cx (Im[i?]) 2 carrying the restoring elastic potential energy of the vacuum 
which acquires an induced shear elasticity in the presence of the charge q [5a,g,h]. 

For stationary state, for which eo(x, t) by definition does not change with time, 



eo(x) — e\tp(x, t)\ 2 = e\tp(x)\ 2 , e= I eo(x)dx = e \tp(x,t)\ 2 dx, 

Jo Jo 



(3) 



where e = L v e E 2 for tp being normalized in [0, L], 
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\tP{x,t)\ 2 dx = 1 (4) 



and thus e — e. The corresponding complex wave function is tp(x,t) = tp{x)e~ l ^ t+a<> \ 
where to = jO, with 7 = 1/ \J\ — w 2 /c 2 [5a,c,d]; ao is an initial phase which affects the field 
configurations of the individual component waves (see further after eqn. 9) but not the total 
particle state (since |e M °| 2 = 1) and can be here simply set to ao = 0. The total linear 
momentum is p = p (x)dx, with po(^) = ^zl^MMMI ■ m the case of V{x) — % being 
constant, p = (e — %)/c. 
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The e, in terms of the total wave ip of a frequency cu, is simultaneously subject to the 
Planck energy equation, 

£ q .n q = n q {b /2n)io, n q = 1, 2, b = h; (for n q = 1) : s q = hw, e = e q ; (5) 

(5) can be obtained as a solution to the energy equation (internal work), of a form analogous 
to the Schrodinger equation, for the harmonically oscillating charge q being point like and yet 
extensive at scale ~ 10~ 18 m, and of a dynamical mass (9H 9 )[5c,j]. Furthermore, for the wave 
train travelling rectilinearly at a finite velocity c which resembles a rigid object [5a, c,c] of a 
finite inertial mass (m), we may employ Newton's law of inertia and obtain po(x) in relation to 
density of mass, mo (a;), (setting f = for defining mass here) and in turn the integration as 

p (x) £ (x) £ 2 £ f L 2 £ p 

m (x) = = — 2— = -z\ip(x,t)\ , m=^ \ip(x,t)\ dx =-$ = -. (6) 

c c z cr cr J cr c 

If ip(x,t) were an arbitrary function of t, the proportionality relation Eq{x) oc \ip(x,t)\ 2 of 
(3) would not generally hold, and s would be generally dependent on t. However under the 
condition of stationary state we assumed, (3) holds generally, as may be readily verified by 
substituting ip{x,t) = tp(x)e^ l ^^ h ^ t in the wave equation sip = Hip of (2), and multiplying it 
by ip* from the left. Its left side gives (3a), or rewritten as 

pL pL 

e (x) = tp*etp = e\tp(x,t)\ 2 (a), e= / ip*iipdx= / tp*Htpdx (b) (7) 

Jo Jo 

where (7b) is given by integrating (7a) on both sides. 

For the specific case of y(x) — and v = 0, (2) has the solution ip = -l= e l ( Kx - nt ) j where 

K = lim„2 / c 2^ k(= lu/c) = Mc/h, J? = lim„2 / c 2^ u = Mc 2 /h and M — lim„2 / c 2^ m[5a,c-e]; 
with this ip in (2) we obtain hflip = h ^ ip. Subtracting the resulting equation from (2) of a 
generally finite V and v, taking limit at v 2 /c 2 — > and at scale of the de Broglie wavevector 
k d = (^)k (<< fc), ^ = lim v 2 /c 2^. ,fe d /K«i i> = &{x)e- mdt , with fi d = we obtain a 

Schrodinger form of wave equation governing the particle's kinetic motion [5b, c], S v ^ = H v ^, 

where S v = i - HQ, H v = H - = + V, and P v = |V. The corresponding energy 

density and integration follow to be 

g vQ (x) =&*i v <P = g v \V(x,t)\ 2 , § v = \ ^*i v ^dx= <P*H v >Pdx, (8) 

Jo Jo 

where § v = S v for \W\ 2 dx = 1; for the case of f{x) = % =constant, S v = \Mv 2 + %. 
We shall proceed the remainder of the discussions in terms of the total wave tp and energy e, 
etc., i.e., the "relativistic" forms, and in one dimension. The conclusions can all be carried in 

2 

a straightforward way over for <P, S v , etc. by taking the limit at ^ — > and scale kd, and to 
three dimensions [5c]. 

Unambiguously, the total energy or mass of the IED particle is as any real particle an 
absolute measure of the existence of the particle. The IED particle of a spatially extensive 
energy and mass as of (3), or (7), and (6) therefore is extensive, here apparently across [0,L] 
(though practically across L v ). And its presence at position x in [0,L] is measured by the 
energy density £o(x), or alternatively mass density mo (a;), which for stationary state contains 
each the common factor \ip\ 2 , denoting by p t (x) as 

Pl (x) = \ip(x,t)\ 2 , Xp, =ip*(x,t)J?ip{x,t) = \ip{x,t)\ 2 J(r. (9) 

where in (9b) p t (x) is written in the more general form of being acted on by a dynamical 
variable operator K. 
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The foregoing descriptions relevant to p l (x) may be summarized as: (i) p t (x) is a distributed 
function in [0,L] describing electromagnetic fields which are implicitly random because the 
initial phase ao of ip (see after eqn. 4) is in general the result of random interactions (Sec. 3). 
(ii) p l (x) is a measure of the fraction of the total energy or mass of the particle presenting at x. 
And £ and m are the p l (x)-weighted integrations across [0,L], and therefore the expectations 
of total energy and mass, (hi) < p,(x) < 1 according to (3) and (4); in particular, p^x) = 
if no portion of the particle is found at x, and p,(x) = 1 if the entire particle is found at x, 
hence a corpuscle. 

Points (i)-(hi) furnish p,(x) = \ip(x,t)\ 2 , or p l (x)=\<P(x,t)\ 2 at scale kd, with all of the 
essential qualifications for p t (x) to be defined as the probability, or probability density here, 
for finding a corresponding portion p l (x) of the IED particle at position xina given stationary 
state. This conclusion is in complete accord with the (mathematical) probability assumption 
in the formal theory of quantum mechanics. 

3. Time evolution of single particle state. Probability density in dynamical 
variable space 

Let now the N = 1 IED particle in the one-dimensional box of volume V — L be subject 
to interactions with its environment, the wall particles and/or environmental radiation fields, 
maintained at a constant temperature T. The particle is here described in the quantum- 
mechanical terms by the wave equation (2) for all regimes, and the description will reduce to 
the classical mechanics limit when the particle's energy levels e v s are on the scale h densely 
nested with one another. Assume that the interactions are weak, infrequent, and elapse over 
a brief time <5i„ only on each occurrence. So the IED particle will maintain a given stationary 
state v (specified by e.g. a total energy, linear momentum, and spin, etc.) for a finite time 
interval t v ,» 5t v , but, subjected to interactions with environment, will explore all possible 
states on a time scale t » t v . 

This may for example describe the realistic case of an electron moving in a fixed orbit about 
the nucleus of an isolated atom until an emission/absorption of radiation, or a conduction 
electron moving at a constant velocity in a small metal specimen until next collision, or a 
spin-half atom traversing a magnetic field with a fixed spin orientation as in the Stern-Gerlach 
experiment (Sec. 3) until externally perturbed. Under the condition of t v » 8t v , contributions 
from transient states may be neglected, and the thermodynamic properties are essentially 
determined by the stationary states. 

Instead of making observations over time, we now "compact" a large H number of the 
states evolved over time (t') into H replicas which each describe the one particle system here 
under condition of constant T, V, N(— 1), and arc uniformly distributed over a total H number 
of possible quantum states of the particle. We thus obtain a canonical ensemble. Each vth 
(eigen) state is associated with an (eig en) wave function ?/v(x,t) — ^v(&)e ^ eu l energy c Vl 
etc., given as the solution of (2), with v = 1, 2, . . . , N. Each vth state occupies a volume element 
in the phase space spanned here by x,p x ; the trajectory of state in this space corresponds to 
the time evolution of state of the actual particle. With (2) being a linear equation, the linear 
sum of tp v (x)s, 



must again be a solution of (2), where c v {t) is a "weight coefficient" or "amplitude" of state v 
and is in general complex. 

That ip ena (x,t) is a linear sum of ip v s implies only the linearity of the wave equation. 
ip cns (x,t) does not describe the total displacement at x due to the superposition of tp v s, since 
the ijj v s occur in reality at different times and never meet. 




(10) 
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Multiplying ip*, (x) on both sides of (10) from the left and integrating over (0, L), we obtain, 
for the ip v s being orthonormal (f^ ip v iip v dx = 5 V ' , Vl 8 v i, v being the Kronecker delta), 

/ ipl(x)ip ens (x,t)dx = / ip*(x)y^c v >(t)ip v >(x)dx = c v {t). (11) 
Jo Jo v , 

For the canonical ensemble of a dimension H assumed large, its average energy, proportional to 
the integral below, is evidently a constant. And the integral itself must aslo be a constant, and 
is normalised to one as 

/ r e n S (x,t)i>en S (x,t)dx = I, Or V C* V C V / \ij> v (*) | 2 dx = V |c„ (t) | 2 = 1 , (12) 

where |c„(i)| 2 = c*(t)c v (t). With for tp in (2), multiplying on both sides of the uth equation 
by c v from the left for v = 1, . . . , H, summing the equations over all u, we obtain a corresponding 
wave equation, 

fVens = ^en S , i = iffy, H = -(fr 2 /m)V 2 + r(z). (13) 

With ?/v for V, making similar algebraic operations as in leading to e of (7), we obtain the 
energy of vth state, ip*(x, Vjih 9 ^"^'^ dx = e v . From the significance of c„ as the amplitude 
of state v in -0 C ns defined in (10), the energy expectation over the states of the ensemble follows 

to be (e v ) = Jq L V'cnsC^i t)ih d ^ on ^ x ^ dx, and 

(£„)=/ y] c* l/ ,(t)ipi,(x)y^ j s u c v (t)'tp u (x)dx = y^s v p v , Pv = \ Cv {t)\ 2 , (u) 

"'° i/' = l i/=l v=l 

where for the final expression for (e v ), \ip^\ 2 dx = 1 is used. The ensemble average (s v ) of (14) 
is equivalent to the time average provided that the IED particle here is ergodic, an assumption 
usually made for real particles. By a similar reasoning as for p l (x) earlier (as p lv (x) = \il> v (x, t)\ 2 
here), the p v — |c„(i)| 2 given in (14) has all of the qualifications to be defined as the probability 
of finding the particle at state v. The function p v is related to e v through a Boltzmann factor, 
oc e~ e ' J l kBT which is a proper problem of statistical mechanics. 

Given an arbitrary dynamical variable K v (x) of an operator J(f v {x) (which may also be an 
independent variable like the x, or the a later) of the single particle for ^th state at position 
x, from the preceding discussions its expectation in both the position and dynamical-variable 
spaces follow to be 

-L oo oo oo „L 

(jf) = c*y>* j(f v {x) c vl ^ V 'dx — ic^i 2 ^"^, ^ v =\ ip*(x)jfciip v (x)dx. (15) 

•*° v=0 v'=0 v=0 ^° 

The tp v s have been introduced in the above as describing events taking place in reality in a 
random fashion one after another in time, only one at a time regardless of being probed by an 
observer or not. This feature is in direct accord with the eigen-state solutions of Schrodinger 
equation for a corresponding system and, as can be expounded case by case, with overall 
experiments. 

A most clear-cut experimental indication is provided by the Stern-Gerlach experiment (Z 
Phys 8, 110, 1921). Here, a silver atom, of a net spin i due to its outermost-shell electron 
and being ordinarily random oriented, is let travel at a constant velocity (v) in the x direction 
across a region d applied with an inhomogeneous transverse (z-direction) magnetic field B z of a 
finite gradient ^§£-{> 0). In the B z field its spin magnetic moment p z = —^g s fi B is quantized 
in orientation (Zeeman effect), a = +1 or — 1. If the atom deterministically maintains a fixed 
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spin orientation, say a = +1, provided also that no other perturbation presents, the eigen state 
function being thus i[> = i[> + , then across the entire d region the atom will constantly be acted 
by a force F z = ^5sMb^^ in the fixed +z direction for a duration t — d/v, and at the exit of 
d be displaced by a finite distance -\-Z\ = ^-^-t 2 in the +z direction (in the experiment the 
atom also traverses an additional free flight distance D which is not of direct concern here). 
Or if a = — 1, the state function being thus ip — tjj-, then F z and z\ change each to the — z 
direction. The above is as observed in the experiment. 

On the contrary, if during each one trip the atom were in a mixture of both spin up and 
down states, so tp = ^ip+ + \^-, then across d the atom would be acted by forces of alternating 
directions. The so produced net displacement at the exit must on average be zero, z\ = 0, 
which is in contradiction with the experimental observation. 



4. The state vector space 



What is mainly new in the foregoing is the concretisation of the physics underlining the 
probabilities in both the position and dynamical variable spaces instrumented by the IED 
model and solutions. The involved mathematical structure and operations are of no essential 
difference from the usual quantum mechanics, since the governing equations at scale kd are the 
common Schrodinger equation. This feature holds through the remainder of discussions. 

As in the unified mathematical formulation of quantum mechanics, the description in the 
(position) coordinate space can be transformed to one in an H-dimensional linear state vector 
space (the Hilbert space for H being infinite) , which for the single IED particle here is spanned by 
the ip„s as coordinate axes, or basic vectors, of a total number H, assuming the tp v s orthonormal 
in (0,L). For many (IED) particles, we have the Fock space. The linear combination of the H 
basic vectors ip„s, V'cns = '52 v Cv'4>v of (10), represents in this space a "(general) state" vector 
which describes based on Sec. 3 the state of the (canonical) ensemble; c v is a component of the 
vector V'ens- And \c v \ 2 is the probability of "finding" the physical system, the single particle in 
(0, L) here, in vth. (stationary) state at a time point t' v on the time evolution axis t' . 

With the ip u s being discrete, the linear mathematical operations in the space spanned by 
the \p v s may be naturally carried out in terms of matrix analysis. If for example ,J^ ai , 
(i, j — 1, 2, . . . , H) are given as two dynamical variables, their two general states are given by 
V'cns. a(x,t) — J2 i ai(t)tpi(x), ip cns .b(x,t) — ^2jbj(t)i^j(x). Their scalar product for example is 
defined by 



(V>a,V>&) 



/ ^a* l/ ipl^2 l b v >ip vl dx = ^alb v 

1/ 7 /' M 



a 2-> • • • > a 'k) 



( h \ 

b 2 



V & H / 



(a\b), (16) 



where J Q ipl^^dx = 5 vu i is used; (tp a — (a\ corresponds in the final expression to a H x 1 matrix 
with its H elements being imaginary and arranged in one raw, and tpb) = \b) a 1 x H matrix 
with its H elements being real and in one column. If a = b = c, (16) returns the result of (12). 
The operator H or ,J(f in general is an H x H matrix. 



5. N weakly interacting particles in a mean interaction field 

Let now N weakly interacting IED particles of the dynamical masses m ai rrib, ■ ■ ., wave functions 
i>a{xi,t), 4>b(x2, t), . . ., and attached with the coordinates xi,X2,- ■ ■, be enclosed in the box of 
size L in an effective mean potential field i / (x 1 ,x 2 , ■ ■ ■ , x N ),= N ^_^ Sill Y^j=i,(j=tt) 
with %j a pairwise interaction potential. Each %jj a {xi,t) is governed by a wave equation given 
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after (2), 

i a il> a (xi,t) = H a tjj a (xi,t), i a = ihd t , H a = -(h 2 /m a )\7 2 + f{xi, . . . , x N ), (17) 

where a = a,b, . . . ,a N , i — 1,2, ... ,N. The N vp a -wa,ves are each distributed across the 
entire L, thereby overlapping with one another. The total field at position Xi is ipj:(xi,t) = 
^2 a ip a {xi,t), which generally contains a finite portion of each of the N waves, or, particles. 
From a measurement of ipz(xi, t) at Xi, we are thus not able to distinguish which particle this is; 
that is, the N particles are indistinguishable, and in E. Schrodinger's terminology (Proc Cam 
Phil Soc 31, 555, 1935), they are entangled with one another. 

All of the expressions up to (16) can accordingly be written down for each ath particle. Of 
specific interest here is the probability density of finding particle a at x%, p la {xi) — \ip a {xi, t)\ 2 
given after (9), with a = a, . . . , ajv, % = 1, . . .,N. The notion of \ip a (xi, t)\ 2 as a probability 
may be advanced for the many-particle system through excising it in a few representative cases 
below, as we shall see, each returning the expected result of the formal quantum mechanics and 
(or) experiment. 

(i) The probability of the simultaneous occurrence of N independent events that particle a 
is at Xi, b is at x 2 , ■ ■ ■ and a N is at x N with the probabilities p l (xi),p t .(£2), ■ • • and p lN (x aN ), 
with the particles being in the mean field effectively non-interacting, is according to the usual 
multiplication rule for probabilities given by 

{a N ,N} 

11 PlA X i) = \^a(.Xi)%l) b {x 2 ) ■ ■ .1p aN (x N )\ 2 . 
{ a ,i} = {a,l} 

Because of indistinguishablity, making simultaneous pair-wise coordinate permutations say from 
x\ to x 2 for particle a and x 2 to x\ for b, leaves the product 1^0(^2)^6(^1)1 unchanged except for 
a change in sign if the permutations are antisymmetric; there are a total AM number of possible 
permutations. To ensure an always positive probability we should introduce as in usual practice 
a "Sign" into the product of N wave functions as (Sign)i/'a(a;i)V'fc( :z: 2) ■ • • 4>a N (x N ). Further 
dividing y/N\ for normalisation, the complete equation for the probability is therefore 

N\ N\ {a N ,N} 

p lN (xi,x 2 , . . . ,x N ) = IVv(xi,x 2 , . • -,x N ;t)\ 2 = II Pi«( x i)> ( 18 ) 

i,V {a,i} = {a,l} 

where 

4> N (x 1 ,x 2 ,...,x N ;t) = (Sign)-j=jp a (x 1 ,t)ip b (x 2 ,t) . . .ip aN (x N ,t), (19) 

with Sign= +1 for t/j N symmetric and —1 antisymmetric in respect to simultaneous permutations 
between a pair of coordinates X{ and X{* . 

One common application of the function form ip N , a complex iV-particle wave function, is 
the description of the simultaneous states of many electrons in an atom, or in the atoms of a 
condensed matter. Multiplying (17) by ^ ^ from the left, summing the equations (17) over 
all as, reorganising, we obtain a wave equation for tp N as 

(iV.ajv) ^ 

e N ^ N = H N i> N , e N = iffy, H N = - ^ V 2 + NY. (20) 

(i,a)=(l,) ^ 

(ii) A beam of (N) incoherent identical particles, each being a travelling wave tjj(xi,t), 
is incident from left onto an opaque screen with two slits A, B, and emerges from these as 
ip A (ri,t),ip B (r 2 ,t), which are then recorded at a point P on a photographic plate on the right. 
We register the passings of the H = 2 slits A and B by the incident beam that have subsequently 
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arrived at P as two events, which are here mutually exclusive and have the probabilities 
P lA (P) — \tp A (P,t)\ 2 and p lB (P) = \ip B (P,t)\ 2 . The probability of the occurrence of at least 
one of the H = 2 two mutually exclusive events is according to the usual addition rule for 
probabilities given by 

PlAB (p) = \MP,t)\ 2 + \MP,t)\ 2 . (21) 

(iii) If the incident particle beam of (ii) consists of (N) coherent waves, then the H = 2 
events are no longer mutually exclusive, but are in a situation as stated by the usual conditional 
probability (p AB = ^(p B p A / B + PaPb/a))', the expressions for p A / B etc. are subject to how the 
events specifically depend on one another. Concretely, for the present case we can readily write 
down the total wave displacement at P as ip AB (P,t) — ip A (P,t) + ip B (P,t). But \ip AB (P,t)\ 2 is 
according to Sec. 2 just the probability p lAB {P), so 

PlAB (p) = |^ AB (P,t)| 2 - \MP,t) + MP,t)\ 2 - (22) 

This will predict the interference fringes observed in double slit experiments. For N identical 
particles as in (ii)-(iii) above, m a = m for all as; directly summing the (17)s over all as gives 
a total wave equation for as ih 9 ^^ 1 '^ — Hiip^(xi,t), Hi = — f^-Vf + V. 

The author expresses thanks to Scientist P.-I. Johansson for his continued moral and private 
funding support for the research, and to a community of national and international distinguished 
physicists for giving their invaluable moral support for the research. 
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